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Abstract 


This paper compares recent spatial anomaly time series of OLR (Outgoing 
Longwave Radiation) and OLRclr (Clear Sky OLR) as determined using CERES and 
AIRS observations over the time period September 2002 through June 2010. We find 
excellent agreement in OLR anomaly time series of both data sets in almost every 
detail, down to the 1x1° spatial grid point level. This extremely close agreement of 
OLR anomaly time series derived from observations by two different instruments implies 
that both sets of results must be highly stable. This agreement also validates to some 
extent the anomaly time series of the AIRS derived products used in the computation of 
the AIRS OLR product. The paper then examines anomaly time series of AIRS derived 
products over the extended time period September 2002 through April 2011. We show 
that OLR anomalies during this period are closely in phase with those of an El Nino 
index, and that the recent global and tropical mean decreases in OLR and OLRclr are a 
result of a transition from an El Nino condition at the beginning of the data record to La 
Nina conditions toward the end of the data period. We show that the relationship 
between global mean, and especially tropical mean, OLR anomalies to the El Nino 
index can be explained by temporal changes of the distribution of mid-tropospheric 
water vapor and cloud cover in two spatial regions that are in direct response to El 
Nino/La Nina activity which occurs outside these spatial regions. 

1. Introduction 

OLR (Outgoing Longwave Radiation) is a critical component of the Earth’s 
radiation budget and represents the total radiation going to space emitted by the earth- 
atmosphere system and integrated over all angles. OLR products have been generated 
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and monitored globally since 1975 based on broad spectral band measurements taken 
at a given satellite zenith angle by the ERB instrument on the Nimbus-6 and Nimbus-7 
satellites (Jacobowitz et al. 1984, Kyle et al. 1993); the ERBE instrument on NOAA-9 
and NOAA-IO; ERBS (Barkstrom 1989), the AVHRR instrument on NOAA operational 
satellites (Gruber et al., 1994 and references therein); and most recently by CERES 
which has flown on EOS Terra since 2000 and EOS Aqua since 2002 (Wielicki et al. 
1996). 

OLR has been widely used as a proxy for tropical convective activity and rainfall, 
particularly in diagnosing and understanding tropical intraseasonal to interannual 
variability and monsoons (e.g., Chiodi and Harrison 2010, Kidson et al. 2002, Hoyos 
and Webster 2007, Barlow et al. 2005, Kiladis et al. 2005, Jones et al. 2004, Wong et al. 
2008, Loeb et al. 2011). In addition, OLR has been used in studies of earth’s radiation 
balance (e.g. Clement and Soden, 2005; Fasullo and Trenberth, 2008) and atmospheric 
model validation (e.g. Allan et al. 2003). More importantly, anomalies and trends of OLR 
have been used to study climate feedbacks and processes (e.g. Chung et al. 2010, 
Huang and Ramaswamy 2009, Soden and Held 2006, Chu and Wang, 1997, Soden et 
al. 2008, Dessler et al. 2008, Dessler 2010, Trenberth et al. 2010). 

OLR at a given location is affected primarily by the earth’s skin surface 
temperature, T S ki n ; skin surface spectral emissivity, e v ; atmospheric vertical temperature 
profile, T(p), and water vapor profile q(p); and the heights, amounts, and spectral 
emissivities of multiple layers of cloud cover. OLR also depends on the vertical 
distribution of trace gases such as 0 3 , CH 4 , C0 2 , and CO. OLR can be computed for a 
specific Field of Regard (FOR), given all the needed geophysical parameters, using an 
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OLR Radiative Transfer Algorithm (RTA). Mehta and Susskind developed such an OLR 
RTA used in conjunction with the TOVS (TIROS Operational Vertical Sounder) retrieval 
methodology (Susskind et al. 1997) in order to generate the TOVS Pathfinder Path-A 
OLR data set (Mehta and Susskind 1999a, 1999b). AIRS OLR is computed using 
AIRS/AMSU sounding products in a completely analogous manner, including use of the 
same Mehta and Susskind OLR RTA (Susskind et al. 2003). 

AIRS measures IR channel radiances over the interval 650 cm' 1 to 2668 cm" 1 . 
Most AIRS results shown in this paper were derived using the AIRS Science Team 
Version-5 retrieval algorithm (Susskind et al. 2011) which generates values of T S kin, e v , 
T(p), q(p), 0 3 (p), cloud parameters, OLR, and OLR C lr valid for each AIRS 45 kmx45 km 
FOR. The AIRS Version-5 OLR product, referred to as F below, is computed as a sum 
of 14 spectral components according to 

F = E}=i Fj = ££(1 - a£ 1} - as 2j ) F jCLS + a£ tj F jiCLD1 + as 2J Fj iCLD2 (1 ) 

where Fj,clr is the computed clear sky flux going to space integrated over all angles 
emanating from spectral band j ; Fj,cwk is the analogous flux emanating from an opaque 
cloud at cloud top pressure pk\ and askj is the radiatively effective cloud fraction which is 
the product of the geometric fractional cloud cover otk for the cloud at pressure pk as 
seen from above and the emissivity of that cloud in spectral interval j. 

Mehta and Susskind parameterize Fj for a given sounding as a function of the 
retrieved surface skin temperature, surface spectral emissivity q in spectral band j, and 
atmospheric temperature, water vapor, and ozone profiles. The parameterization 
coefficients used by Mehta and Susskind are computed based on line-by-line 
calculations (Susskind and Searl, 1978) which use the atmospheric line parameter data 
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base of McClatchey etal. (1972). The spectral intervals used in Equation 1 range from 
2 cm' 1 through 2750 cm" 1 . There is no need to make radiometric measurements at all 
frequencies in order to perform the calculation shown in Equation 1. The surface 
spectral emissivity s v is determined as a function of frequency over the AIRS spectral 
range using AIRS observations. Surface emissivities at frequencies lower than 650 cm' 1 
are set equal to those at 650 cm' 1 and are irrelevant with regard to the computation of 
OLR in any event because the atmosphere is opaque at those frequencies. The AIRS 
Version-5 retrieval algorithm determines the effective cloud fraction askj at 800 cm" 1 for 
each of two cloud layers k. ask and is assumed to be independent of frequency in the 
calculation of OLR. No other approximations are made in the calculation of Equation 1. 
AIRS OLRclr, the clear sky OLR, is also a computed product obtained using Equation 1 
but setting both asi and as 2 equal to zero. Geophysical parameters are determined from 
AIRS observations under both cloud-free and cloudy conditions, though their quality is 
poorer under very cloudy conditions. The CERES OLRclr product provides OLR for 
CERES footprints observed to be cloud-free according to coincident MODIS spectral 
radiance measurements. The MODIS cloud mask used by CERES is described in 
Minnis et al. (201 1 ). 

2. Overview of the Paper 

This paper has two main thrusts. The first thrust of the paper is its comparison of 
values of OLR and OLRclr determined using CERES and AIRS observations, and their 
anomaly time series, over the time period September 2002 through June 2010. This 
comparison shows that the CERES and AIRS OLR anomaly time series are in very 
close agreement with each other. Agreement of OLR anomaly time series, obtained by 
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two different instruments in two very different manners, is a very significant finding as it 
lends credence to the scientific validity of the results obtained from each instrument. A 
close agreement of OLR anomaly time series derived using two different approaches 
does not necessarily imply agreement of the OLRclr anomaly time series, however, 
because OLRclr is more difficult to obtain via either approach. 

With regard to AIRS, OLRclr is a computed product using the retrieved 
geophysical parameters. Retrieved geophysical parameters are generally of poorer 
accuracy under very cloudy conditions, especially at or near the surface. For this 
reason, the AIRS Version-5 OLRclr product is computed only for cases in which the 
AIRS retrieved cloud fraction is less than 90% and which also pass an OLRclr quality 
control procedure which indicates the retrieval is of acceptable accuracy down to the 
surface (Susskind et al. 2011). We produce successful AIRS Version-5 OLRclr 
products in roughly 75% of the scenes observed by AIRS. OLR is always computed 
from AIRS data, both because of the need for complete sampling, and also because 
computed values of OLR are not affected significantly by surface and atmospheric 
conditions beneath the cloud in very cloudy cases. With regard to CERES, spurious 
values of OLRclr could occur if a scene which is considered to be clear actually 
contained some residual cloud cover. The largest factor that would negatively affect the 
comparison of OLRclr anomaly time series obtained from AIRS and CERES results 
from the sampling differences between the two sets of cases included in each 
ensemble. 

Both AIRS and CERES OLR data sets show that a linear fit through their 
respective global mean monthly mean anomaly time series over the period of 
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September 2002 through April 2010 contains a substantial negative slope, on the order 
of -0.075 W/m 2 /yr, and an even larger decrease of tropical mean OLR on the order of 
-0.1 W/m 2 /yr. There are very large spatial variations of the changes in OLR in the 
tropics, with local values ranging frorr2.8 W/m 2 /yr to +3.1 W/m 2 /yr. We observe a 
correlation of 0.95 between the spatial patterns of changes found in the in AIRS and 
CERES OLR data sets during this time period. In addition, slopes of the linear least 
squares fits of monthly mean anomaly time series averaged over different spatial 
regions agree on the order of ±0.01 W/m 2 /yr and the spatial-temporal correlation of the 
two sets of the OLR anomaly time series in the tropics is 0.993. This extremely close 
agreement of OLR anomaly time series derived from observations by two different 
instruments, and determined in totally independent and different manners, implies that 
both sets of results must be highly stable. This level of agreement strengthens the 
overall confidence in the ability of spaceborne sensors to detect large scale global 
signals and allowing the creation of more reliable multi-sensor, merged data sets. This 
agreement also validates to some extent the AIRS derived products used in the 
computation of the AIRS OLR product. Furthermore this agreement indicates we can 
use the anomaly time series of AIRS derived products to explain the factors contributing 
to the anomaly time series of both the AIRS and CERES OLR products. 

The second thrust of the paper examines anomaly time series of AIRS derived 
products over the extended time period September 2002 through April 2011. We show 
OLR anomalies during this period are closely in phase with those of an El Nino index, 
and the global and tropical mean decreases in OLR and OLRclr are a result of a 
transition from an El Nino condition at the beginning of the data record to La Nina 
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conditions at the end of the data period. The relationship between global mean, and 
especially tropical mean, OLR anomalies to the El Nino index can be explained by 
temporal changes in two spatial regions of the distribution of mid-tropospheric water 
vapor and cloud cover that are in direct response to El Nino/La Nina activity which 
occurs outside of these regions. One of these regions is east of the area of El Nino/La 
Nina activity, and the second region is south-southwest of that area. Changes in 
otherwise global mean and tropical mean OLR, but computed in areas excluding these 
two regions, are very close to zero over the time period under study. 

3. AIRS and CERES OLR data sets used 

In this paper we primarily use the operational monthly mean OLR and OLRclr 
data products derived using the AIRS and CERES Science Team methodologies 
respectively. We obtained the AIRS OLR from the Goddard DISC and the CERES 
products from the Langley ASDC. AIRS was launched on the EOS Aqua satellite in a 
1:30 AM/PM local crossing time orbit in May 2002. The operational processing of AIRS 
data began after AIRS became stable in September 2002. We use the AIRS Version-5 
monthly mean Level-3 1°x1° latitude-longitude grid products which contain separate 
products generated for each of the 1 :30 AM and PM local time orbits. We averaged the 
AM and PM products together to generate and use a single monthly mean product on 
the 1° x 1° grid for each month included in the data set. In addition to AIRS OLR and 
OLRclr, we also use AIRS Level-3 T S kin, q(p), and as cloud products to analyze the 
factors contributing to the anomaly time series of OLR and OLRclr. 

CERES has flown on both EOS Terra, which was launched in December 1999 on 
a 10:30 AM/PM local crossing time orbit, as well as EOS Aqua, the same platform that 
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carries AIRS. The CERES Science Team generates a number of different OLR data 
sets using CERES observations. The latest versions of the longest record CERES OLR 
data sets are referred to as the CERES SSFIdeg-lite Edition-2.5 data sets, which like 
AIRS, are presented on 1°x 1° latitude-longitude grid. At the time of this writing, the 
CERES Terra Edition-2.5 OLR data sets extended to June 2010, the CERES Aqua 
Edition-2.5 OLR data sets extended to February 2010, and the AIRS Level-3 products 
extended to April 2011. CERES SSFIdeg-lite Edition-2.5 uses the latest calibration 
improvements (Priestley et al. 2011) with Edition-2 CERES cloud retrievals (Minnis et al. 
2008, Minnis et al. 2011), angular dependence models (Loeb et al. 2005), and time- 
space averaging (Young et al. 1998). 

This paper also shows some results comparing OLR computed using the 
prototype AIRS Version-6 Science Team retrieval algorithm with the Version-5 OLR 
product computed for the same days. AIRS Version-6 uses an improved OLR RTA 
(lacono et al. 2008) in the computation of OLR. This new RTA has two important 
upgrades compared to Mehta and Susskind (1999a,b). Most significantly, the new OLR 
RTA is generated using more up to date line absorption parameters, especially for the 
very strong water vapor absorption band near 300 cm' 1 . In addition, the new OLR 
calculation allows for inclusion of the effects of variations in CO 2 concentration over 
time, as well as those of other minor absorption species such as CO, CH4, and N2O, in 
the calculation of OLR. The Version-5 OLR RTA did not include these effects and 
parameterized atmospheric transmittances only in terms of variable temperature, water 
vapor, and ozone. The AIRS Version-6 retrieval algorithm also has other improvements 
in methodology which lead to improved values of the geophysical parameters 
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themselves. The AIRS Science Team Version-6 retrieval algorithm is expected to 
become operational in late 201 1 . We obtained the OLR results shown in the paper from 
JPL. They are not available to the public at this time. 

4. Comparison of AIRS and CERES OLR and OLRclr Data Records 

Figure la shows global monthly mean values of AIRS OLR, as well as CERES 
Aqua and CERES Terra Edition-2.5 OLR (henceforth referred to as CERES Aqua and 
CERES Terra OLR) for the period starting September 2002 and extending until the end 
of each of the current data sets. This figure, as well as subsequent figures, contains no 
smoothing unless otherwise noted. Figure 1b shows analogous results for OLRclr. The 
AIRS OLR time series contains the symbol * on the months November 2003 and 
January 2010. AIRS OLR and OLRclr products for parts of these months were missing 
from the daily AIRS data record. We generated AIRS monthly mean OLR products for 
these months synthetically, on a grid box basis, by setting grid point differences 
between AIRS and CERES OLR for an incomplete month equal to the average value of 
the corresponding AIRS/CERES difference for the previous and subsequent month. We 
used the same procedure to generate the AIRS OLRclr data records for those two 
months. We use these synthetic OLR and OLRclr monthly mean records as if they 
were actual observations in the subsequent discussion. 

Some monthly mean CERES OLRclr data points for individual 1° x 1° grid boxes 
have also been synthesized because either no CERES OLRclr monthly mean data 
existed for these grid points or the data values were outliers when compared to AIRS 
OLRclr values. For each grid box for each month, we first eliminated any CERES 
OLRclr value that differed by more than 20 W/m 2 from the corresponding AIRS OLRclr 
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value. Roughly 2% of the CERES OLRclr monthly mean gridded values were 
eliminated using this criterion. We then spatially interpolated the remaining values of the 
difference, CERES OLRclr minus AIRS OLRclr, to generate synthetic values for the 
missing grid boxes. These synthetic differences were then added to the AIRS OLRclr 
product for the corresponding grid box to generate the synthetic value of CERES 
OLRclr which is used in all subsequent calculations. 

We observe a number of features apparent from Figures la and 1b. The most 
prominent result is that to first order, the AIRS and CERES OLR data sets appear to be 
biased compared to each other, as are the AIRS and CERES OLRclr data sets. Figures 
2a and 2b show the differences between the monthly global mean values of OLR and 
OLRclr shown in Figures la and 1b, as well as the difference between CERES Terra 
and CERES Aqua. Figures 2a and 2b also contain dashed lines showing the average 
value of each difference. AIRS Version-5 OLR shows a nearly constant bias, with an 
average value of 9.05 W/m 2 , compared to CERES Terra and 9.49 W/m 2 compared to 
CERES Aqua OLR. We observe a small essentially repetitive seasonal cycle, with 
maxima and minima in June and December respectively, in the difference between 
AIRS and CERES Terra OLR. Part of this annual cycle in the difference between AIRS 
and CERES OLR may be the result of the large diurnal cycle of OLR over land. CERES 
tries to capture the diurnal cycle by fitting the data at the Terra and Aqua overpass 
times (twice per day at low latitudes) with a half-sine fit (Young et al. 1988). AIRS 
averages daytime and nighttime OLR observations together but does not make any 
other correction for diurnal cycle. The difference between AIRS and CERES Aqua OLR 
is similar to that between AIRS and CERES Terra OLR, but somewhat more variable 
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over time. CERES Terra OLR is slightly larger than CERES Aqua OLR, especially 
before January 2005. As with AIRS OLR, there is a small mostly repetitive annual cycle 
in the difference between CERES Terra OLR and CERES Aqua OLR from 2005 on. The 
differences between AIRS and CERES OLRclr are similar to, but smaller than, those of 
OLR, with regard to their mean value and their seasonal cycle. The mean value of AIRS 
minus CERES Terra OLRclr is 6.37 W/m 2 , which is roughly 2.7 W/m 2 less than that of 
AIRS minus CERES Terra OLR. 

The large biases between AIRS and CERES OLR and OLRclr data records are 
at first disconcerting but are readily understood. The AIRS OLR product derived using 
the AIRS Science Team Version-6 retrieval algorithm, which will become operational in 
late 2011, is expected to have much smaller biases compared to CERES Terra OLR 
than does AIRS Version-5 OLR. The nearly constant bias between OLR as computed 
based on AIRS products and observed by CERES is primarily a result of the use of an 
old set of line by line absorption coefficients in the parameterization of the Version-5 
OLR RTA (Mehta and Susskind 1999a,b), compared to the improved OLR RTA (lacono 
et al. 2008) used in Version-6. The main difference between the two OLR 
parameterizations is that lacono et al. (2008) has more absorption in the water vapor 
rotational band than does Mehta and Susskind (1999a). This change would tend to 
generate lower values of OLR, especially under very moist conditions. 

We tested the latest version of the AIRS Version-6 retrieval algorithm on seven 
days in different years and different seasons, ranging from September 6, 2002 through 
May 20, 2010. The difference in OLR between the two versions has significant spatial 
variability even on a given day. In most cloudy areas, which are less sensitive to 
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absorption by tropospheric water vapor, Version-6 OLR tends to be higher than 
Version-5 OLR. Conversely Version-6 OLR tends to be considerably lower than 
Version-5 OLR in most regions away from clouds, especially for moister cases. Figure 
3a shows the spatial distribution of the difference between the 7-day average of OLR 
computed using Version-6 and computed using Version-5. Figure 3b shows analogous 
results for OLR C lr. The most important point to notice is that the Version-6 global mean 
OLR averaged over seven days is lower by 7.06 W/m 2 than the Version-5 global mean 
OLR product. This indicates that AIRS Version-6 global mean OLR should agree with 
CERES Edition-2.5 OLR to within about 2 W/m 2 . Figure 3b shows that global mean 
Version-6 OLRclr for these seven days is 7.54 W/m 2 less than Version-5 OLRclr. This 
suggests that AIRS Version-6 OLRclr will also be very close to CERES OLRclr on the 
average. 

Table la shows values of the slopes of the linear least squares fits of the 
differences among AIRS, CERES Terra and CERES Aqua OLR, as well as the mean 
differences and standard deviations of these OLR time series over the period 
September 2002 through August 2009. Table 1b shows analogous results comparing 
the OLRclr time series. The comparison statistics are given for this portion of the total 
overlap time period because it represents differences over a complete 7-year period. It 
is important that these slopes be compared over complete years to minimize artifacts 
due to effects of the annual cycle differences in each data set. 

The slope of the differences between AIRS and CERES OLR and OLRclr 
records over this 7-year time period is not affected to first order by a constant bias and 
gives a preliminary indication of the stability of different data sets with regard to each 
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other over the seven year period under study. The slope for the time series AIRS minus 
CERES Terra OLR is considerably smaller than that of AIRS minus CERES Aqua OLR. 
The standard deviation between the AIRS OLR and CERES Terra OLR time series is 
also considerably smaller than that between AIRS OLR and CERES Aqua OLR, and 
closer to that between CERES Terra OLR and CERES Aqua OLR. These same relative 
results also hold when comparing time series of OLR C lr. This does not necessarily 
imply that the CERES Terra OLR record is better than the CERES Aqua OLR record. 
Nevertheless, based on the results shown in Tables la and 1b, and the fact that the 
CERES Terra OLR data set extended further in time than the CERES Aqua OLR data 
set, we chose to use the CERES Terra OLR and OLRclr products in the OLR 
comparisons between CERES and AIRS shown in the remainder of this paper. 

5. Comparison of AIRS and CERES Terra OLR Anomaly Time Series 

In this section, we compare anomaly time series of AIRS and CERES OLR and 
OLRclr and their Average Rates of Change (ARC's) computed over different domains: 
global mean; tropical mean; zonal mean; and grid point by grid point. Anomalies 
represent differences of monthly mean values of a product from their climatological 
value as determined by each instrument. Therefore, biases between instrumental 
records do not affect the relationship between their anomalies to first order. We define 
the ARC for a given data set as the slope of the linear least squares fit passing through 
an anomaly time series. We use the term Average Rate of Change to describe the 
slopes of anomaly time series rather than the term Trend, which is generally used to 
characterize long-term multi-decadal data sets rather than the seven plus year period 
studied in this paper. The results show that even though there is a significant bias 
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between AIRS and CERES data records, there is still an excellent agreement between 
details of AIRS and CERES OLR and OLRclr anomaly time series in almost all details. 

We generated the monthly mean OLR and OLRclr climatologies on a 1°x 1° 
spatial resolution for each month by taking the average of the grid box value for that 
month over a 7-year time period. OLR and OLRclr anomalies for a given month in a 
given year, on a 1° x 1° spatial grid, are defined as the difference between their monthly 
mean values and their monthly climatologies for that grid box. The area mean anomaly 
for a given month is the cosine latitude weighted average of the grid box anomalies 
contained in the area under consideration (Global, Tropical, North America etc.). 

5.1 . Global Mean and Tropical Mean Anomaly Time Series 

This section compares global mean and tropical mean anomaly time series of 
AIRS and CERES OLR and OLRclr- Changes in global mean quantities over time are 
often used as an indication of whether the earth-atmosphere system has been changing 
as a function of time (e.g., the discussion of "global warming"). Such a discussion can 
be misleading, both because changes in the earth-atmosphere system over time are not 
spatially homogeneous, and also because changes of a parameter over a limited time 
period are in no way indicative of what may occur over later time periods. We also 
examine tropical mean anomaly time series. The tropics, defined as the band 20°N 
through 20°S, constitutes 34% of the area of the globe. As such, features found in 
tropical monthly mean anomalies can contribute significantly to those found in global 
monthly mean anomalies. 

We show that global and tropical mean anomalies of AIRS and CERES OLR are 
highly correlated in time with each other, as well as with the AIRS El Nino index. The 
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AIRS El Nino index represents the monthly mean oceanic sea surface temperature 
(SST) anomaly determined using the AIRS Version-5 T S kin product, averaged over the 
spatial area 15°N to 15°S latitude and 140°W westward to 160°E longitude. The reason 
for the selection of this El Nino area will be discussed later in the paper. We also show 
reasonably good agreement between AIRS and CERES anomaly time series of 
OLRclr, both with each other and with the El Nino index. The latter is an important 
finding because it demonstrates that while changes in cloud cover in response to El 
Nino/La Nina activity are playing an important role in the recent decreases of global 
mean, and especially tropical mean OLR, other short-term changes influenced by El 
Nino have also taken place that affect OLRclr. 

Figure 4a shows the global anomaly time series of AIRS Version-5 OLR and 
CERES Terra OLR for the periods September 2002 through April 201 1 and September 
2002 through June 2010, respectively, as well as the difference between the two sets of 
monthly mean anomalies for the overlap time period. Figure 4b shows analogous 
results for tropical mean OLR anomalies, and Figures 5a and 5b show analogous global 
mean and tropical mean anomaly time series for the AIRS and CERES OLRclr 
products. Figures 4b and 5b also include the El Nino index, multiplied by 3 in Figure 4b 
and by 2 in Figure 5b. Tropical mean OLR and OLRclr anomalies tend to track those of 
the El Nino index in phase fairly closely. Amplitudes of the largest tropical mean OLRclr 
anomalies also match those of the El Nino index multiplied by 3, and those of OLRclr 
anomalies match the El Nino index multiplied by 2. This shows that the longest tropical 
mean OLR anomalies are roughly 50% larger than the largest tropical mean OLRclr 
anomalies. Positive values of the El Nino index (2003, 2005, 2007, early 2010) 
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correspond to El Nino periods, and negative values (2008, late 2010) correspond to La 
Nina periods. 

Table 2a shows the global and tropical mean values of the Average Rate of 
Change for AIRS OLR and CERES Terra OLR anomalies over the time period 
September 2002 through June 2010, the standard deviations between the two sets of 
global and tropical anomaly time series, and the temporal correlations between each 
global and each tropical anomaly time series. The uncertainties shown in this table and 
all other tables represent the 99% confidence interval. The agreement of the Average 
Rates of Change of both global mean and tropical mean anomaly time series found in 
the AIRS and CERES OLR records is very good, with a difference of about 
0.02 W/m 2 /yr. The standard deviation of the difference between the two global OLR 
anomaly time series is considerably smaller than that between the two OLR time series 
because the effects of the small essentially constant annual cycle in the differences 
between AIRS and CERES Terra OLR have been removed to first order in the 
generation of the anomaly time series. In addition to considerations relating to the 
diurnal cycle of OLR, a part of the small annual cycle in the difference between AIRS 
and CERES OLR can be attributed to some extent to the fact that there is more water 
vapor in the atmosphere, and hence a larger error in the Version-5 OLR calculations, in 
the Northern Hemisphere summer than in the Northern Hemisphere winter. The 
temporal correlations of the two global and tropical OLR anomaly time series are 0.953 
and 0.983 respectively. Both show that global mean OLR has decreased on the 
average on the order of -0.075 W/m 2 /yr over the common time period September 2002 
through June 2010, and tropical mean OLR has decreased at a rate of roughly 
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-0.1 W/m 2 /yr from the beginning of the time period to the end. The close agreement of 
global and tropical mean ARC’S of AIRS and CERES OLR anomaly time series is more 
significant than the values themselves. 

Figures 4a and 4b show that an onset of negative mean OLR anomalies began in 
late 2007, with tropical mean values generally considerably larger than global mean 
anomalies, especially after mid-2007. The decrease in global mean OLR in late 2007 is 
strongly influenced by the significant reduction in tropical mean OLR which started a few 
months earlier. Tropical mean OLR anomalies, and to a lesser extent global mean OLR 
anomalies, became positive starting in late 2009, roughly coincident with the onset of 
another El Nino event. The AIRS OLR product, which extends to April 2011, shows that 
very substantial negative global and tropical mean OLR anomalies occurred in the 
period starting mid-2010, in which a substantial La Nina event occurred which continues 
through the time period observed by AIRS shown in this paper. As a result of inclusion 
of these large negative OLR anomalies, the negative Average Rate of Change of global 
mean AIRS OLR over the extended time period through April 201 1 increased to -0.093 
W/m 2 /yr, compared to -0.088 W/m 2 /yr through June 2010, and the negative ARC of the 
tropical mean AIRS OLR increased to -0.189 W/m 2 /yr compared to -0.1 13 W/m 2 /yr. 

Figure 4b shows that the anomaly time series of tropical mean OLR closely 
follows that of the El Nino index. The ARC of the El Nino index, computed over the 
extended time period through April 2011, is -0.11 K/yr, which indicates that on the 
average, the El Nino region has been cooling over the last 8 years. The global mean 
and tropical mean OLR anomaly time series indeed correlate very highly with each 
other as well as with the El Nino index. 
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Table 2c shows temporal correlations between global mean and tropical mean 
anomaly time series of OLR and OLRclr as well as correlations of anomaly time series 
with the AIRS El Nino index. Correlations using AIRS time series are shown above the 
diagonal in bold and those using CERES time series are shown beneath the diagonal. 
The temporal correlation of the CERES global and tropical mean OLR anomaly time 
series is 0.55, and the corresponding correlation for the AIRS time series is 0.56. This 
shows that tropical anomalies provide a significant contribution to the global OLR 
anomaly time series. The CERES and AIRS tropical mean OLR anomaly time series 
also correlate highly with the El Nino index, with temporal correlations of 0.68 and 0.69, 
respectively. Both sets of global OLR anomaly time series also show moderate 
correlations with the El Nino index which are smaller than those found in the tropics. 
These correlations of global and tropical anomaly time series with the El Nino index 
imply that the recent short term decreases in global and tropical OLR are the result of 
changes from El Nino conditions at the beginning of the time series to La Nina 
conditions at the end. 

One might not expect as good agreement between anomaly time series of 
CERES OLRclr as with OLR, if for no other reason than there are significant sampling 
differences between the cases included in each OLRclr data set. Given this caveat, it is 
worthwhile to compare both OLRclr anomaly time series to see the extent that features 
found in one data set corroborate the same features found in the other. It is also 
important to see the extent that features in the OLRclr anomaly time series are 
consistent with those of OLR. A considerable agreement between OLR and OLRclr 
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anomalies indicates that a substantial part of the OLR anomalies can be attributed to 
something other than clouds. 

Table 2b shows that the AIRS and CERES OLRclr anomaly time series are 
highly correlated with each other, though the correlation is not as high as between the 
OLR time series. Part of the decrease in correlation is a result of the reduction in values 
of each set of OLRclr anomalies compared to those of OLR. In addition, the standard 
deviations of the OLRclr anomaly differences are somewhat larger than those of OLR. 
Perhaps more significant is the fact that AIRS and CERES global mean and tropical 
mean OLRclr anomaly time series are each very highly correlated with those of OLR as 
shown in Table 2c. Moreover, each tropical OLRclr anomaly time series has roughly 
comparable correlations with the El Nino index as does tropical OLR. This indicates that 
a considerable portion of global and tropical OLR anomalies are influenced by El 
Nino/La Nina activity and arise from something other than changes in cloud cover. 

While AIRS and CERES OLRclr anomaly time series are highly correlated with 
each other, Figures 5a and 5b both show a positive drift in time of the OLRclr anomaly 
differences. Table 2b shows that the global and tropical mean values of the ARC’S of 
AIRS OLRclr are very close to zero, while those of CERES OLRclr are very close to 
those of CERES OLR. This apparent discrepancy may well be a result of the significant 
sampling differences found in the two different OLRclr data sets. 

5.2. The Spatial Distribution of ARC’S of OLR and OLRci r 

This section compares the spatial distribution of ARC’S of AIRS and CERES OLR 
and OLRclr with each other. First, we compare zonal mean ARC’S, and then compare 
ARC’S on a°1latitude by 1° longitude basis. These comparisons not only show 
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excellent agreement of ARC’S of AIRS and CERES OLR products on a small spatial 
scale, but also depict the spatial regions that have been contributing to the short term 
decreases in global mean and tropical mean values of OLR over the period under study. 

Figure 6a shows the zonal mean ARC’S of the AIRS and CERES Terra OLR 
records for the period September 2002 through June 2010, as well as the difference of 
the two sets of zonal mean ARC’S. There is again excellent agreement between the 
zonal mean structures of ARC’S in both OLR data sets, which have a latitudinal 
correlation of 0.97. It is apparent from Figure 6a that the majority of the decrease in 
global OLR during the time period under study originates in the tropics south of 8°N. 
Other areas of negative OLR ARC’S occur near 60°S and 60°N latitudes. OLR 
increased over this time period north of 70°N. AIRS negative zonal mean OLR ARC’S 
are somewhat larger than those of CERES between 5°S and 45°S and AIRS OLR 
ARC’S are more positive (or less negative) than those of CERES poleward of 70°. 

Figure 6b shows analogous results for OLR C lr. With regard to AIRS, the 
structure of zonal mean ARC’S of OLR and OLRclr are generally similar to each other, 
with a correlation of 0.86. The agreement is especially good poleward of 60° latitude. 
AIRS zonal mean OLRclr ARC’S in the tropics are significantly smaller than those of 
OLR. In addition, there is a major qualitative difference between OLR and OLRclr near 
5°N latitude, which has a very large negative ARC of AIRS OLR and essentially a zero 
ARC of AIRS OLRclr- ARC’S of CERES zonal mean OLRclr are very close to those of 
AIRS OLRclr poleward of 60°N and 60°S latitudes, but are generally on the order of 0.1 
W/m 2 /yr more negative between 6CfS and 50°N. The latitudinal correlation of AIRS and 
CERES ARC’S of zonal mean OLRclr is 0.92, which is slightly lower than that found for 
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AIRS and CERES OLR. The latitudinal correlation of ARC’S of CERES OLR and 


OLRclr is 0.85, which is again similar to the correlation of zonal mean ARC’S of AIRS 
OLR and OLRclr. The fact that zonal mean ARC’S of OLR and OLRclr are highly 
correlated with each other shows that much of the cause of the zonal mean structure of 
OLR ARC’S is coming from something in addition to changes in cloud cover. 

The spatial distributions of global OLR ARC’S over the time period September 
2002 through June 2010 are shown in Figures 7a and 7b for AIRS and CERES, 
respectively. Both fields are presented on the 1° latitude by 1° longitude grid on which 
the data sets are given. The precise values of the ARC’S shown in Figures 7a and 7b 
are not necessarily meaningful, both because of the limited statistical significance of the 
value for each grid point taken over a seven year ten month time period, and also 
because as discussed previously, these values are dependent on the time period over 
which they are calculated. More significant than the values of the ARC’S shown in 
Figure 7 is their very coherent spatial structure and its implication with regard to 
atmospheric processes. 

Figures 7a and 7b demonstrate two very important points. The first is the virtually 
indistinguishable spatial distributions of the ARC’S of AIRS OLR and of CERES OLR. 
Figure 7c shows this difference, with a spatial correlation of 0.948, and a standard 
deviation of 0.18W/m 2 /yr. The global mean AIRS OLR ARC for this period is 
0.023 W/m 2 /yr lower (more negative) than that of CERES Terra. This small difference is 
not monolithic, but occurs primarily near 30°S latitude, especially over Eastern Australia 
and South Africa, both of which contain large negative OLR ARC’S. 
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The most important point of Figures 7a and 7b is that not only do OLR Average 
Rates of Change contain a pronounced zonal mean structure, there is considerable 
longitudinal structure at given latitudes as well. This longitudinal structure is particularly 
noteworthy in the tropics. As shown in Figure 7, even though zonal mean tropical OLR 
ARC’S are very negative, indicative of increased zonal mean precipitation over the time 
period September 2002 through June 2010, positive OLR ARC’S as large as 
3.2 W/m 2 /yr, indicative of a mean decrease in precipitation, exist in the vicinity of the 
equatorial dateline. These are more than compensated for, in the zonal mean sense, by 
negative OLR ARC’S at other longitudes, as large a§.9 W/m 2 /yr near the equator 
over Indonesia in the vicinity of 120°E longitude. 

Figures 7a-c, and some subsequent figures, contain rectangles surrounding the 
area between 8°N to 20°S and 140°W eastward to 10°E, which we will refer to as OLR 
Region 1, and the area between 15°S to 30°S and 140°E eastward to 160°W, which we 
will refer to as OLR Region 2. ARC’S of OLR within both these areas are very negative 
over the period September 2002 through June 2010. The decreases in OLR that took 
place in these two regions during the period under study will be shown later in the paper 
to be the major cause of the negative values of the global and tropical mean OLR 
ARC’S shown in Table 1. OLR Regions 1 and 2 were selected entirely based on the 
results shown in Figure 7 and were chosen both so as to be rectangular and also to 
encompass key features shown in both Figures 7a and 7b. OLR Region 2 is essentially 
in the “heart” of the area referred to as the South Pacific Convergence Zone, which 
varies its location according to the phase of the El Nino - Southern Oscillation (ENSO) 
and of the Interdecadal Pacific Oscillation (Folland et al. 2002; Brown et al. 201 1 ). 
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Figures 8a-c show the spatial distribution of ARC’S of AIRS and CERES OLRclr 
over the time period September 2002 through June 2010, as well as their difference. 
Regions 1 and 2 are also indicated in these figures. The general agreement of the 
spatial distribution of ARC’S of AIRS and CERES OLRclr is again very good, with a 
spatial correlation of 0.762. Figure 8c shows that much of the positive bias of AIRS 
zonal mean ARC’S of OLRclr compared to CERES between 50°N and 60°N latitude 
occurs over the oceans. At high latitudes, both the patterns and magnitudes of ARC’S of 
OLR and OLRclr agree well. In the tropics, however, the patterns of ARC’S of OLR and 
OLRclr agree well, but the magnitudes of ARC’S of OLR are considerably larger than 
those of OLRclr- This indicates that changes in cloud cover over the time period under 
study are playing a major role with regard to changes in tropical OLR, but are playing a 
much smaller role with regard to changes in extra-tropical OLR. 

5.3. Longitudinal Distribution of Equatorial Anomaly Time Series: Hovmoller 
Diagrams 

Figures 7 and 8 show that the tropics contain large spatially coherent areas with 
alternating values of positive and negative ARC’S of OLR and OLRclr over the time 
period under study. This section compares Hovmoller diagrams of the longitudinal 
distribution of tropical anomaly time series of CERES and AIRS OLR and OLRclr- 
Agreement of AIRS and CERES tropical Hovmoller diagrams is significant because it 
demonstrates that there is not only agreement in ARC’S of OLR and OLRclr, but in the 
tropical OLR anomaly time series themselves. These anomaly time series play a very 
important role in explaining the tropical distribution of the ARC’S of OLR and OLRclr 
shown in Figures 7 and 8. 
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Figures 9a and 9b present Hovmoller diagrams, showing time series of monthly 
mean OLR anomalies (vertical scale), integrated over the latitude range 5°N through 
5°S, in each 1° longitude bin (horizontal scale) for the periods for which AIRS and 
CERES OLR records exist: September 2002 - April 201 1 for AIRS and September 2002 
- June 2010 for CERES. The difference between these two figures in the overlap time 
period is shown in Figure 9c. Figures 9a-9c, and all subsequent Hovmoller diagrams, 
have a small amount of smoothing applied to them. A five point (5 month) smoothing 
was applied in the vertical and a fifteen point (15 degree) smoothing was applied in the 
horizontal to minimize the effects of small discontinuities between adjacent rectangular 
grid points on the figures. Most of the region covered is ocean. There are three 
relatively small land areas near the equator: South America, Africa, and Indonesia. 
These land areas each lie between the three sets of narrow vertical lines shown in 
Figure 9. There are also two other vertical lines at 140°W longitude and 10°E longitude, 
encompassing the longitudinal range contained within Region 1. 

The two sets of Hovmoller diagrams are essentially identical, with a correlation 
coefficient of 0.993 between them. Some of the largest differences between the AIRS 
and CERES tropical anomaly time series occur in November 2003 and January 2010, 
the two months for which AIRS data was synthesized. These differences would have 
been much larger if the AIRS “monthly mean” OLR products stored at the Goddard 
DISC were used in the calculations, because in both cases the AIRS “monthly mean” 
products represented averages over less than a month time period while the CERES 
data represented observations taken over the entire month. Indeed, when the AIRS 
data contained at the DISC for these two months was used originally in the generation 


24 



of Figure 9, those two months showed pronounced signals in the original Figure 9c. This 
alerted us to check, and correct for, the cause of this problem. 

The anomaly time series shown in Figures 9a and 9b depict the phase 
relationship of OLR anomalies in different longitudes as a function of time. Such figures 
in turn provide insight into the spatial distribution of tropical ARC’S in the vicinity of the 
equator such as shown in Figures 7 and 8. In the longitudinal band 16(yw westward to 
140°E, equatorial OLR anomalies were very negative in late 2002/early 2003, 
corresponding to a period of increased precipitation, and were very positive from mid- 
2007 through early 2009 corresponding to decreased precipitation. This gave rise to the 
substantial positive OLR ARC shown in Figure 7 over the regionIS - 5°S, 160°W to 
140°E. Figure 9a shows that AIRS OLR anomalies in this same region are also very 
positive starting mid-2010. Therefore, OLR ARC’S over the extended time period 
September 2002 - April 2011 would be even more positive in this area than those 
shown in Figure 7a covering a shorter time period. Figure 7 shows very negative values 
of OLR ARC’S near the equator between 100E and 140°E longitudes. Figure 9 shows 
that equatorial OLR anomalies between 100E and 140°E are out of phase with those 
between 140°E and 160°W and are of comparable magnitude. Figures 9a and 9b also 
show that equatorial OLR anomalies 16tW eastward to 40°E , within the longitudinal 
domain of OLR Region 1, tend to be smaller than, and out of phase with, those from 
160°W westward to 140°E. This gives rise to the negative equatorial OLR ARC’S shown 
in Figure 7 in OLR Region 1 . 

Figures lOa-c show analogous results for OLR C lr. The scale of Figure 10 is one 
third that of Figure 9, so as to better highlight the smaller OLRclr anomalies. Hovmoller 


25 



diagrams of AIRS and CERES OLRclr anomaly are again in excellent agreement, with 
a correlation of 0.913 between them. The patterns of OLRclr anomalies are also very 
similar to those of OLR, albeit with smaller values of the anomalies, which are roughly 
one-third the value of the OLR anomalies. Hovmoller diagrams of AIRS OLR and 
OLRclr have a spatial-temporal correlation of 0.793 with each other over the time 
period September 2002 - June 2010, and those of CERES OLR and OLRclr have a 
correlation of 0.697. This again demonstrates that while temporal changes in cloud 
cover are playing a very large role in tropical OLR anomalies, other factors which are in 
phase with changes in cloud cover play a significant role as well. 

6. The Effect of Phases of El Nino/La Nina on Tropical Water Vapor, Cloud 
Cover, and OLR Anomaly Time Series 

Figures 4 to 10 show that the spatial patterns of the AIRS and CERES Average 
Rates of Change of OLR and OLRclr over the time period September 2002 through 
June 2010 are in excellent agreement with each other, as are their anomaly time series 
averaged over different spatial regions. Both CERES and AIRS OLR products show that 
the common time period under study is marked by a substantial decrease in global 
OLR, on the order of0.075 W/m 2 /yr, averaged over the globe and). 10 W/m 2 /yr 
averaged over the tropics. This agreement of Average Rates of Change of OLR 
anomaly time series derived from observations by two different instruments, in totally 
independent and different manners, implies that both sets of OLR products must be 
stable over the 7 year 10 month period in which they were compared. There should be 
little question that there actually was a decrease of Global mean OLR on the order of 
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0.075 W/m 2 /yr over the time period September 2002 through June 2010, and that the 
majority of the decrease occurred in the tropics. 

This result, found by both CERES and AIRS, should not be taken as indicative in 
any way as to what will happen in the future. It mainly shows that OLR anomalies and 
their Average Rates of Change can be determined very accurately by two totally 
independent instrumental and theoretical approaches. The agreement of anomaly time 
series of OLR as observed by CERES and computed from AIRS derived products also 
indirectly validates the anomaly time series of the AIRS derived products used in the 
computation of AIRS OLR. Moreover, it further indicates that anomaly time series of 
AIRS derived products can be used to explain the factors contributing to anomaly time 
series of OLR. This section of the paper uses the anomaly time series of the AIRS 
derived products and their ARC’S over the extended time period September 2002 - April 
2011 to explain the factors contributing to the anomaly time series of OLR and OLRclr 
over that time period. 

Figure 11a shows the spatial distribution of the Average Rate of Change of the 
AIRS Version 5 surface skin temperature over the period September 2002 - April 201 1 . 
A number of important features are found in Figure 1 la. While the global mean surface 
skin temperature ARC is essentially zero over this time period, there are areas where 
significant positive and negative T S kin ARC’S exist. Both polar regions have been 
warming over this time period, especially near the North Pole. There is considerable 
warming and cooling structure in areas over Northern Hemisphere extra-tropical land, 
and there has been substantial cooling over much of Africa, especially south of 15°S, as 
well as over much of Australia. The tropics are marked by a substantial oceanic cooling 
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over the region 15°N to 15°S and 140°W longitude westward to 160°E which we will 
refer to as the AIRS El Nino region. The AIRS El Nino region, indicated by the green 
rectangle in Figure 11a, was defined to encompass the region of large negative sea 
surface temperature ARC’S shown in the figure. There has also been substantial 
oceanic warming over this time period in the areas to the south, west, and north of the 
AIRS El Nino region. The AIRS El Nino region as defined in this paper is similar to the 
commonly used Nino 4 region, 5°N - 5°S, 120°W westward to 170°W (Yeh, et al. 2009). 

Figure 11b shows the global distribution of ARC’S of 500 mb specific humidity 
(%/yr) over the time period September 2002 to April 2011. OLR is very sensitive to the 
concentration of mid-upper tropospheric water vapor in very moist (i.e. tropical) areas, in 
the sense that increasing water vapor concentration increases atmospheric absorption 
in some spectral regions and therefore lowers OLR, everything else being equal. The 
area of strong equatorial surface temperature cooling between 160°E and 140°W, 
shown in Figure 11a, called the AIRS El Nino region, is marked by a significant mid- 
tropospheric drying during this period, and the surrounding areas of warming surface 
skin temperature are marked by considerable mid-tropospheric moistening during the 
period under study. There is also considerable mid-tropospheric moistening from 8°N to 
20°S in the Atlantic Ocean and especially in the Eastern Pacific Ocean over this time 
period in areas where there has been at most a very small increase in surface skin 
temperature. In all of these regions, the changes in atmospheric water vapor, rather 
than in surface skin temperature, have been more significant contributors to the ARC’S 
of AIRS OLR and OLR C lr over the extended time period September 2002 - April 201 1 , 
shown in Figures 11c and lid. Figures 1 1 b-1 1 d all contain the boundaries of OLR 
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Regions 1 and 2 defined previously. OLR Regions 1 and 2 both correspond to locations 
in which there has been considerable mid-tropospheric moistening over this time period, 
as well as considerable decreases in both OLR and OLRclr. Figures 11c and lid are 
similar to, but different than, Figures 7a and 8a which show analogous ARC’S computed 
over the shorter time period extending only to June 2010. As stated previously, the 
AIRS global mean ARC’S of OLR and OLRclr over the extended time period are -0.093 
W/m 2 /yr and -0.018 W/m 2 /yr respectively. Both values are considerably more negative 
than they were over the shorter time period as a result of inclusion of the latest La Nina 
time period during which there were substantial negative global mean anomalies of OLR 
and OLRclr. 

Figures 11c and lid show that the spatial distribution of ARC’S of AIRS OLR and 
OLRclr over extra-tropical land areas closely track those of T S kin with areas of warming 
skin temperatures corresponding to areas of increasing OLR and OLRclr and vice 
versa for cooling areas. This demonstrates the expected result that a warming surface 
should result in increasing flux to space, all other things being held constant. 

The spatial relationship between ARC’S of OLR and T S kin in the tropics is quite 
different than in the extra-tropics, especially over ocean. There has been a significant 
equatorial cooling between September 2002 and April 2011 in the El Nino region 
surrounding the equator from 160°E eastward to 140°W. This is a reflection of the fact 
that the time period under study started with an El Nino event in 2002, and had La Nina 
events in 2007 and 2010. This change from El Nino conditions to La Nina conditions is 
the driving force behind the recent negative global mean and tropical mean OLR ARC’S. 
The AIRS El Nino region of T S kin cooling is marked by significant increases in both OLR 
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and OLRclr. This region is surrounded to the north, south, and west by areas of 
oceanic warming. All of these oceanic areas are marked by large negative OLR ARC’S. 
It is extremely important to note that equatorial OLR also shows large negative ARC’S 
over OLR Region 1, 8°N to 20°S, 140°W eastward to 10°E, in which there are no 
spatially coherent changes with regard to T S kin- 

Figure lid shows that spatial ARC’S of clear sky OLR, which to first order are 
independent of changes in cloud cover, closely match ARC’S of OLR poleward of EjO 
both in phase and magnitude. This shows that OLR ARC’S at high latitudes are not 
influenced substantially by changes in cloud cover ae. In the tropics, the patterns of the 
ARC’S of OLR and clear sky OLR are also similar to each other, but the magnitudes of 
the relative ARC’S are substantially larger with regard to OLR than OLRclr. Tropical 
OLRclr ARC’S are dominated by those of 500 mb specific humidity q 5 oo- Tropical OLR 
ARC’S are also affected by changes in water vapor, but are dominated primarily by 
those of cloud cover ae. ARC’S of fractional cloud cover (not shown) closely match 
those of 500 mb specific humidity. It is the combination of the effects of anomalies and 
trends of both tropospheric water vapor and fractional cloud cover that give rise to the 
anomalies and ARC’S of OLR as observed by CERES and confirmed by AIRS. 

Table 3 shows the spatial correlation of ARC’S of OLR, OLRclr, ore, q 50 o, and T skin 
over the time period September 2002 through April 201 1 . Spatial correlations of ARC’S 
of these quantities over the latitudinal domain 60N through 90°N are shown above the 
diagonal in bold, and analogous spatial correlations over the latitudinal domain°N5 
through 15°S are shown below the diagonal. At high latitudes, shown in bold, ARC’S of 
OLR and OLRclr have a spatial correlation of 0.83 with each other. Both sets of ARC’S 
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are also highly positively correlated with those of T S km, especially OLRclr. At high 
latitudes, ARC’S of specific humidity and fractional cloud cover are not highly correlated 
with each other or with any of the other geophysical parameters shown in the table. 

Spatial correlations of ARC’S of the same geophysical parameters in the tropics, 
shown below the diagonal, are very different from those at high latitudes. The 
correlations of ARC’S of OLR and OLRclr in the tropics is again very high (0.80). In the 
North Polar region, the high correlation of ARC’S of OLR and OLR clr arose because 
the spatial distribution of ARC’S of OLR and OLRclr were both highly positively 
correlated with those of T S kin- In the tropics, this high correlation of OLR and OLRclr 
arises because both are highly negatively correlated with ARC’S of fractional cloud 
cover ae and 500 mb specific humidity q 50 o- Unlike at high latitudes, ARC’S of tropical 
Tskin are not highly correlated with those of any other of the geophysical parameters 
studied. Tropical correlations of ARC’S and OLR and OLRclr with those of T S kin are 
actually somewhat negative because tropical correlations of T sk i n with those of ae and 
q 50 o are somewhat positive. The high correlation of ARC’S of OLRclr with those of ae, 
on which OLRclr does not depend, arises because ARC’S of ae and q 50 o are 
themselves highly positively correlated with each other. 

Figure 1 2a shows the Hovmoller diagram of monthly mean T S km anomalies for the 
period September 2002 through April 2011. The vertical green lines delineate the 
longitudinal band included in the El Nino region, 160°E longitude eastward to 140°W 
longitude. The largest SST anomalies tend to occur on either side of the dateline 
between these longitudinal limits. Figure 12a demonstrates that the large negative SST 
ARC near the equator extending from about 160°E eastward to 140°W is the result of 
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the transition from an El Nino at the end of 2002 to La Ninas over the time periods late 
2007 through 2008, and especially late 2010 through early 2011. Equatorial T S kin 
anomalies between 100°E and 140°E tend to be smaller than, and of opposite sign to, 
those in the vicinity of the dateline. This gives rise to the band of weaker positive SST 
ARC’S near the equator from 100°Eto 140°E. 

Figure 12b shows the Hovmoller diagram of 500 mb specific humidity. The 
vertical gray lines delineate the longitudinal extent of OLR Region 1, extending from 
140°W eastward to 10°E. 500 mb specific humidity anomalies near the dateline in 
general follow those of SST very closely both in magnitude and in phase. This results 
from the fact that positive SST anomalies in the El Nino region correspond to periods of 
increased convection in that area, leading to enhanced moisture of the mid-troposphere. 
Conversely, negative SST anomalies in the El Nino area correspond to periods of 
decreased convection (increased subsidence) leading to periods of a drier mid- 
troposphere. Water vapor anomalies over Indonesia, from roughly 100°E to 140°E, are 
not only out of phase with those near the dateline, as are those of T S km, but are also of 
comparable magnitude to water vapor anomalies near the dateline. This is the result of 
the westward shift of the area of maximum convection during La Nina periods from the 
dateline to Indonesia. This out of phase relationship gives rise to the very substantial 
positive ARC of 500 mb specific humidity over Indonesia during this time period, as 
shown in Figure 1 1b. 

Figure 11b also shows substantial positive qsoo ARC’S in some tropical locations 
in which no significant changes in T S km exist. The most notable of these is off the west 
coast of South America, in the vicinity of 8°N to 20°S from 140°W eastward to 80° W, 
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which is a part of OLR Region 1. There is also another region of positive qsoo ARC’S 
near the equator going across South America and extending eastward along the Atlantic 
Ocean to about 10°E longitude, also contained within OLR Region 1. Figure 12b shows 
that equatorial water vapor anomalies off the west coast of South America tend to be 
out of phase with those at the dateline, especially during the large La Nina events in 
2007-2008 and 2010-201 1 . This indicates that La Nina periods of decreased convection 
near the dateline correspond to periods of increased convection eastward of HOE , 
which is the result of the eastward shift of the convective branch of the Walker 
circulation during La Nina periods (Power and Smith 2007; Zhou et al. 2011). The same 
relationship is found to a lesser extent over the Atlantic Ocean extending to ICE 
longitude at the eastern end of OLR Region 1 . 

Figure 12c shows the Hovmoller diagram of the AIRS effective cloud fraction. 
The anomalies of effective cloud fraction averaged over 3SI to 5°S latitude are closely 
in phase with those of 500 mb specific humidity, with a temporal/spatial correlation of 
0.717. This correlation is a result of the fact that periods of increased convection (or 
subsidence) correspond to increases (or decreases) of not only mid-tropospheric water 
vapor but also cloud cover, especially with regard to high clouds. Anomalies of OLRclr 
are driven primarily by those of qsoo, and those of OLR are driven primarily by effective 
cloud fraction, are, with a smaller contribution from q 50 o- The high correlation between 
tropical anomalies of OLR and OLRclr are the result of the correlation between 
anomalies of are and q 50 o- 

Table 4 shows the spatial-temporal correlations of the Hovmoller diagrams of 
anomalies of OLR, OLRclr, ere, q 50 o, and T S kin derived from AIRS observations over the 


33 



period September 2002 through April 2011. Unlike Table 3, this table is symmetric 
above and below the diagonal because only one set of correlation statistics is 
presented. Spatial-temporal anomaly correlations between OLR, OLR C lr, ere, q 5 oo, and 
Tskin averaged over the latitude range °5I through 5°S are extremely close to those 
shown in the lower half of Table 3, representing spatial correlations of ARC’S of these 
geophysical parameters in the tropics. This shows that not only are the spatial patterns 
of Average Rates of Change of tropical OLR, OLR C lr, are, and q 50 o highly correlated with 
each other, but the detailed structure of their spatial anomalies as a function of time are 
also highly correlated. 

7. Attribution of Recent Decreases in Global and Tropical Mean OLR to 
Changes Contained within OLR Regions 1 and 2 

Figure 1 1c shows that the largest OLR ARC’S occur in the tropics over a region 
covering Indonesia on the one hand, and near the dateline on the other. OLR ARC’S in 
these areas are roughly equal to each other and of opposite sign, as are the OLR 
anomalies shown in Figure 9a. The effects of the large positive and negative tropical 
OLR ARC’S near the dateline and over Indonesia tend to cancel in the zonal mean 
sense. As shown in Figure 6a, the negative tropical mean ARC comes primarily from 
the latitudinal range 8°N to 20°S. This in turn originates primarily from the subset of the 
tropics contained within Region 1, 8°N through 20°S and 140°W eastward to 10°E, 
outlined in Figures llb-d. Analogous results are found in Figure 6b with regard to zonal 
mean ARC’S of OLRclr. 

Figure 13a shows the AIRS OLR anomaly time series averaged over OLR 
Region 1 in dark blue. Figure 13a also includes in black the El Nino index, this time 


34 



multiplied by 4, rather than multiplied by 3 as was done in Figure 4b. The Region 1 OLR 
anomaly time series closely follows that of the El Nino index, but is lagged in time (El 
Nino occurs first) by about 3 months. Therefore, in Figure 13a, Region 1 OLR 
anomalies are plotted 3 months earlier than they actually occurred. For this reason, the 
dark blue line starts 3 months before the black line, and also ends 3 months earlier than 
the black line. The correlation coefficient of the 3 month lagged Region 1 OLR anomaly 
with that of the El Nino index is 0.898. Figure 13a also contains the anomaly time series 
of OLRclr averaged over Region 1 , shown by the light blue line, which like the dark blue 
line is lagged 3 months from the El Nino index. OLRclr anomalies in Region 1 are 
smaller than those of OLR, but are highly correlated with those of OLR, with a 
correlation coefficient of 0.873. The 3 month lagged OLRclr anomalies are also highly 
correlated with those of the El Nino index, with a correlation of 0.821 . 

The OLR and OLRclr anomalies averaged over OLR Region 1 are considerably 
larger than corresponding tropical mean anomalies shown in Figures 4b and 5b. 
Consequently the negative ARC’S of OLR and OLRclr averaged over OLR Region 1, 
which are -0.530 W/m 2 yr and -0.150 W/m 2 /yr respectively, are considerably larger than 
the corresponding negative ARC’S of tropical mean OLR and OLRclr given by -0.188 
W/m 2 /yr and -0.074 W/m 2 /yr. These values are included in Table 5, which shows area 
mean Average Rates of Change of the AIRS anomaly time series of OLR and OLRclr, 
computed over the time period September 2002 through April 2011, averaged over 
different spatial domains. 

The negative tropical mean OLR ARC over the period September 2002 through 
April 2011 computed as previously, but after replacing OLR ARC’S in OLR Region 1 by 
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zeroes, is reduced from -0.188 W/m 2 /yr to -0.032 W/m 2 /yr, and the negative global 
mean OLR ARC is reduced from -0.093 W/m 2 /yr to -0.039 W/m 2 /yr. When an analogous 
procedure is done with regard to OLRclr, the negative tropical mean OLRclr ARC is 
reduced from -0.074 W/m 2 /yr to -0.030 W/m 2 /yr and the negative global mean ARC of 
OLRclr is reduced from -0.018 W/m 2 /yr to -0.003 W/m 2 /yr. This indicates that a large 
part of the recent negative global mean and tropical mean OLR and OLRclr ARC’S 
results from the contribution of OLR and OLRclr anomalies within OLR Region 1 . These 
anomalies are in turn in phase with the El Nino index, but lagged in time by roughly 3 
months. It is for this reason that the global and tropical OLR and OLRclr anomaly time- 
series shown in Figures 4 and Figure 5 are in phase with, and highly correlated with, El 
Nino/La Nina activity as noted earlier. 

Figures 11c and lid also show that there are substantial negative ARC’S of OLR 
and OLRclr within the rectangular box surrounding the area from 15°S to 30°S and 
westward from 150°W to 140°E, referred to as OLR Region 2. One-third of this region is 
in the tropics and the remaining two-thirds is in the subtropics. OLR Region 2 is also 
marked by large positive ARC’S of mid-tropospheric water vapor, as shown in Figure 
1 1 b, and cloud cover (not shown) over the period under study. Figure 1 3b shows in dark 
blue and light blue the time series of OLR anomalies and OLRclr anomalies in OLR 
Region 2, superimposed on the time series of the El Nino index, this time multiplied by 
10, without a time lag. As shown in Table 5, OLR in Region 2 has an extremely large 
negative ARC of -1 .74 W/m 2 /yr, and OLRclr in Region 2 has a very large negative ARC 
of -0.54 W/m 2 /yr. All three curves in Figure 13b are highly correlated with each other in 
time. The correlation between the unlagged anomaly time series of the El Nino index 
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and OLR in region 2 is 0.854, the correlation between the El Nino index and the 
anomaly time series of OLRclr in Region 2 is 0.781, and the correlation between the 
anomaly time series of OLR and OLRclr in region 2 is 0.936. As shown in Table 5, the 
ARC’S of global and tropical mean OLR and OLRclr, computed after zeroing out OLR 
and OLRclr ARC’S in both OLR Region 1 and OLR Region 2, are reduced to essentially 
zero. This shows that the large negative global and tropical ARC’S of OLR and OLRclr 
over the period September 2002 through April 201 1 can be attributed almost completely 
to temporal changes in mid-tropospheric water vapor and cloud cover that have taken 
place primarily within OLR Region 1 along the equator in the eastern Pacific and 
Atlantic Ocean, with a further contribution from OLR Region 2 primarily in the southern 
hemisphere subtropics south-southwest of the El Nino region, covering Eastern 
Australia and parts of the southwestern Pacific Ocean. 
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Table la. Comparison of Global OLR Time Series 
September 2002 through August 2009 



AIRS Minus 
CERES Terra 

AIRS Minus 
CERES Aqua 

CERES Terra Minus 
CERES Aqua 

Slope [W/m 2 /yr] 

-0.0195 

-0.0690 

-0.0495 

Mean [W/m 2 ] 

9.07 

9.53 

0.46 

Standard Dev. [W/m 2 ] 

0.345 

0.475 

0.243 


Statistical comparisons of time series for the period September 2002 through August 
2009 for AIRS minus CERES Terra OLR, AIRS minus CERES Aqua OLR, and CERES 
Terra minus CERES Aqua OLR. Shown are the slopes of the least squares linear fits, 
and mean differences and standard deviations of the different OLR time series. 


Table 1b. Comparison of Global OLRclr Time Series 
September 2002 through August 2009 



AIRS Minus 
CERES Terra 

AIRS Minus 
CERES Aqua 

CERES Terra Minus 
CERES Aqua 

Slope [W/m 2 /yr] 

0.0475 

0.0762 

0.0287 

Mean [W/m 2 ] 

6.34 

6.49 

0.15 

Standard Dev. [W/m 2 ] 

0.325 

0.339 

0.513 


Statistical comparisons of time series for the period September 2002 through August 
2009 for AIRS minus CERES Terra OLRclr, AIRS minus CERES Aqua OLR C lr, and 
CERES Terra minus CERES Aqua OLRclr. Shown are the slopes of the least squares 
linear fits, and mean differences and standard deviations of the different OLR time 
series. 
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Table 2a. OLR Anomaly Time Series Comparison 
September 2002 through June 2010 


Data Set 

Global 

Tropical 

AIRS ARC (W/m^/yr) 

-0.0881 ±0.0055 

-0.1125 ±0.0128 

CERES Terra ARC (W/m 2 /yr) 

-0.0655 ± 0.0050 

-0.1006 ±0.0121 

AIRS Minus CERES STD (W/m 2 ) 

0.147 

0.210 

AIRS/CERES Correlation 

0.953 

0.983 


Global and Tropical statistical comparisons of OLR anomaly time series for the period 
September 2002 through June 2010 for AIRS and CERES Terra OLR. Shown are the 
Average Rates of Change, the standard deviations between the anomaly time series, 
and the temporal correlations of the anomaly time series. 


Table 2b. OLRclr Anomaly Time Series Comparison 
September 2002 through June 2010 


Data Set 

Global 

Tropical 

AIRS ARC (W/m'/yr) 

-0.0045 ± 0.0036 

-0.0177 ±0.0067 

CERES Terra ARC (W/m^/yr) 

-0.0703 ± 0.0043 

-0.1018 ±0.0078 

AIRS Minus CERES STD (W/m 2 ) 

0.195 

0.273 

AIRS/CERES Correlation 

0.857 

0.922 


Global and Tropical statistical comparisons of OLRclr anomaly time series for the 
period September 2002 through June 2010 for AIRS and CERES Terra OLRclr. Shown 
are the Average Rates of Change, the standard deviations between the anomaly time 
series, and the temporal correlations of the anomaly time series. 
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Table 2c. Correlations between OLR and OLRclr Anomaly Time Series in Different 

Domains 

AIRS and CERES 



Global 

OLR 

Tropical 

OLR 

Global 

OLRclr 

Tropical 

OLRclr 

El Nino 
Index 

Global 

OLR 

— 

0.56 

0.77 

0.48 

0.47 

Tropical 

OLR 

0.55 

— 

0.51 

0.87 

0.69 

Global 

OLRclr 

0.79 

0.56 

— 

0.65 

0.46 

Tropical 

OLRclr 

0.55 

0.87 

0.73 

— 

0.63 

El Nino 
Index 

0.41 

0.68 

0.52 

0.70 

— 


Temporal correlations of AIRS and CERES OLR and OLRclr global and tropical 
anomaly time series. Correlations using AIRS data records are shown above the 
diagonal in bold and those using CERES data are shown beneath the diagonal. 


Table 3. Spatial Correlations of ARC’S of AIRS Derived Products over the Time Period 

September 2002 through April 201 1 
90°N to 60°N and 15°N to 15°S 



OLR 

OLRclr 

ae 

9500 

Tskin 

OLR 

— 

0.83 

0.07 

0.23 

0.77 

OLRclr 

0.80 

— 

0.43 

0.36 

0.91 

ae 

-0.91 

-0.80 

— 

0.11 

0.44 

qsoo 

-0.81 

-0.82 

0.69 

— 

0.33 

Tskin 

-0.45 

-0.08 

0.26 

0.43 

— 


Spatial correlation of ARC’S of AIRS derived products over the time period September 
2002 through April 201 1 . Spatial correlation for the region 60 through 90°N are shown 
in bold above the diagonal. Spatial correlations for the region IS through 15°S are 
shown beneath the diagonal. 
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Table 4. Spatial-temporal Correlations of AIRS Anomaly Hovmoller Diagrams 5°N-5°S 
September 2002 Through April 201 1 



OLR 

OLRclr 

ae 

qsoo 

Tskin 

OLR 

1.000 

0.80 

-0.93 

-0.79 

-0.45 

OLRclr 

0.80 

1.000 

-0.77 

-0.72 

-0.05 

ae 

-0.93 

-0.77 

1.000 

072 

0.30 

qsoo 

-0.79 

-0.72 

0.72 

1.000 

0.53 

Tskin 

-0.45 

-0.05 

0.30 

0.53 

1.000 


Spatial-temporal Correlations of AIRS Anomaly Hovmoller Digrams 5 -5°S 
September 2002 Through April 2011. Values above and below the diagonal represent 
the same phenomena and are equal to each other. 


Table 5. Area Mean Average Rates of Change of OLR and OLRclr (W/m 2 /yr) 

September 2002 through April 201 1 


Spatial Area 

OLR 

OLRclr 

Global 

-0.093 ± 0.0045 

-0.018 ±0.0031 

Tropical 

-0.188 ± 0.01 1 1 

-0.074 ± 0.0059 

Region 1 

-0.530 ± 0.0229 

-0.150 ±0.0097 

Region 2 

-1.740 ±0.0672 

-0.540 ±0.0218 

Global excluding Region 1 

-0.039 ± 0.0035 

-0.003 ± 0.0026 

Tropical excluding Region 1 

-0.032 ± 0.0080 

-0.030 ± 0.0038 

Global excluding Region 1 and 2 

-0.004 ± 0.0035 

0.008 ± 0.0026 

Tropical excluding Region 1 and 2 

0.002 ± 0.0078 

-0.018 ±0.0035 


Area mean Average Rates of Change of AIRS OLR and OLRclr anomaly time series 
over the period September 2002 through April 201 1 computed over different spatial 
domains. 
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Figure and Table Captions 
Figure la 

Global monthly mean OLR from Version-5 AIRS and Edition-2.5 CERES aboard Aqua 
and Terra for the period September 2002 through April 2011. 

Figure 1b 

Global monthly mean OLR C lr from Version-5 AIRS and Edition-2.5 CERES aboard 
Aqua and Terra for the period September 2002 through April 2011. 

Figure 2a 

Global monthly mean OLR differences among AIRS, CERES Terra, and CERES Aqua; 
and dashed line showing the average value of each OLR difference for the overlap 
periods between each data set. 

Figure 2b 

Global monthly mean OLRclr differences among AIRS, CERES Terra, and CERES 
Aqua; and dashed line showing the average value of each OLRclr difference for the 
overlap periods between each data set. 

Figure 3a 

The AIRS Science Team Version-6 OLR minus Version-5 OLR for the 1:30 AM/PM 
average of the seven day period, September 6, 2002, January 25, 2003, September 29, 
2004, August 5, 2005, February 24, 2007, August 10, 2007, and May 30, 2010. 
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Figure 3b 

The AIRS Science Team Version-6 OLRclr minus Version-5 OLR clr for the 1:30 
AM/PM average of the same seven days. 

Figure 4a 

Global mean AIRS and CERES Terra monthly mean OLR anomaly time series and their 
anomaly differences for the overlap period September 2002 through June 2010. 

Figure 4b 

Tropical (20°N to 20°S) mean monthly mean OLR anomaly time series and anomaly 
differences of AIRS and CERES Terra for the overlap period September 2002 through 
June 2010. Also shown in black is the AIRS El Nino index multiplied by 3. 

Figure 5a 

Global mean monthly mean OLRclr anomaly time series and anomaly differences of 
AIRS and CERES Terra for the overlap period September 2002 through June 2010. 

Figure 5b 

Tropical (20°N to 20°S) mean monthly mean OLRclr anomaly time series and anomaly 
differences of AIRS and CERES Terra for the overlap period September 2002 through 
June 2010. Also shown in black is the AIRS El Nino index multiplied by 2. 
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Figure 6a 

Zonal mean Average Rates of Change (ARC's) for AIRS and CERES Terra OLR 
monthly mean time series for the period September 2002 through June 2010, as well as 
the difference of the two sets of ARC’S. The majority of the decrease in global OLR 
during this period originates in the tropics south of 8°N. 

Figure 6b 

Zonal mean ARC'S for AIRS and CERES Terra OLR C lr monthly mean time series for 
the period September 2002 through June 2010, as well as the difference of the two sets 
of ARC’S. 

Figure 7a 

Spatial 1° latitude by 1° longitude distribution of global AIRS OLR ARC’S over the time 
period September 2002 through June 2010. OLR Region 1 is outlined in gray and OLR 
Region 2 is outlined in black in this and most subsequent figures showing spatial 
distributions of ARC’S of different parameters. 

Figure 7b 

Spatial 1° latitude by 1° longitude distribution of global CERES OLR ARC’S over the 
time period September 2002 through June 2010. 

Figure 7c 

Spatial distribution of ARC’S of AIRS minus CERES OLR. 
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Figure 8a 

Spatial 1° latitude by 1° longitude distribution of global AIRS OLRclr ARC’S over the 
time period September 2002 through June 2010. 

Figure 8b 

Spatial 1° latitude by 1° longitude distribution of global CERES OLRclr ARC’S over the 
time period September 2002 through June 2010. 

Figure 8c 

Spatial distribution of ARC’S of AIRS minus CERES OLRclr. 

Figure 9a 

Hovmoller diagram for time series of monthly mean AIRS OLR anomalies integrated 
over the latitude range 5°N through 5°S in each 1° longitude bin for the period 
September 2002 through April 2011. The longitudinal domain of OLR Region 1 is 
indicated by the gray vertical lines in this and most subsequent Hovmoller diagrams. 

Figure 9b 

Hovmoller diagram for time series of monthly mean CERES OLR anomalies integrated 
over the latitude range 5°N through 5°S in each 1° longitude bin for the period 
September 2002 through June 2010. 
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Figure 9c 

Hovmoller diagram for the difference between AIRS OLR and CERES OLR. 

Figure 10a 

Hovmoller diagram for time series of monthly mean AIRS OLRclr anomalies integrated 
over the latitude range 5°N through 5°S in each 1° longitude bin for the period 
September 2002 through April 201 1 . 

Figure 10b 

Hovmoller diagram for time series of monthly mean CERES OLRclr anomalies 
integrated over the latitude range 5°N through 5°S in each 1° longitude bin for the 
period September 2002 through June 2010. 

Figure 10c 

Hovmoller diagram for the difference between AIRS OLR and CERES OLRclr 
F igure 11a 

Spatial distribution of Average Rate of Change of AIRS surface skin temperature (K/yr) 
over the period September 2002 through April 201 1 . The spatial domain of the AIRS El 
Nino region is outlined in green. 
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Figure 1 1b 

Spatial distribution of Average Rate of Change of AIRS 500 mb specific humidity (%/yr) 
over the period September 2002 through April 2011. 

Figure 11c 

Spatial distribution of Average Rate of Change of AIRS OLR (W/m 2 /yr) over the period 
September 2002 through April 2011. 

Figure lid 

Spatial distribution of Average Rate of Change of AIRS Clear Sky OLR (W/m 2 /yr) over 
the period September 2002 through April 201 1 . 

Figure 12a 

Hovmoller diagram of AIRS Surface Skin Temperature anomalies (K). The longitudinal 
domain of the El Nino region is enclosed between the green vertical lines. 

Figure 12b 

Hovmoller diagram of AIRS 500 mb specific humidity anomalies (%). The longitudinal 
domain of Region 1 is enclosed between the gray vertical lines. 

Figure 12c 

Hovmoller diagram of AIRS effective cloud fraction anomalies (%). The longitudinal 
domain of Region 1 is enclosed between the gray vertical lines. 
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Figure 13a 

AIRS OLR and OLRclr anomaly time series averaged over OLR Region 1 

superimposed on the AIRS El Nino index multiplied by 4 and plotted 3 months earlier. 

Figure 13b 

AIRS OLR and OLRclr anomaly time series averaged over OLR Region 2 

superimposed on the AIRS El Nino index multiplied by 10 without a time lag. 

Table 1 

Statistical comparisons of anomaly time series for the period September 2002 through 
August 2009 for AIRS minus CERES Terra OLR, AIRS minus CERES Aqua OLR, and 
CERES Terra minus CERES Aqua OLR. Shown are the slopes of the least squares 
linear fits, and mean differences and standard deviations of the different OLR time 
series. 

Table 2 

Global and Tropical statistical comparisons of OLR anomaly time series for the period 
September 2002 through February 2010 for AIRS and CERES Terra OLR. Shown are 
the Average Rates of Change, the standard deviations between the anomaly time 
series, and the temporal correlations of the anomaly time series. 
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Table 3 


Spatial correlation of ARC’S of AIRS derived products over the time period September 
2002 through April 2011. Spatial correlations for the region 60 N through 90N are 
shown in bold above the diagonal. Spatial correlations for the region 15N through 15°S 
are shown beneath the diagonal. 

Table 4 

Spatial-temporal correlations of AIRS anomaly Hovmoller Diagrams 5°N-5°S September 
2002 through April 2011. Values above and below the diagonal represent the same 
phenomena and are equal to each other. 

Table 5 

Area mean Average Rates of Change of AIRS OLR and OLR C lr anomaly time series 
over the period September 2002 through April 2011 computed over different spatial 


domains. 



September 2002 through April 201 1 
Global OLR Time Series (Watts/m 2 ) 



AIRS Version-5 
CERES Terra 
CERES Aqua 

Interpolated AIRS Version-5 

Figure la. Global monthly mean OLR from Version-5 AIRS and Edition-2.5 CERES 
aboard Aqua and Terra for the period September 2002 through April 2011. 
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September 2002 through April 201 1 
Global Clear Skv OLR Time Series (Watts/m 2 ) 



AIRS Version-5 

CERES Terra 
CERES Aqua 

Interpolated AIRS Version-5 


Figure 1b. Global monthly mean OLR C lr from Version-5 AIRS and Edition-2.5 CERES 
aboard Aqua and Terra for the period September 2002 through April 201 1 . 
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September 2002 through June 2010 
Global OLR Time Series Differences (Watts/m 2 ) 



AIRS Version-5 minus CERES Terra Edition 2.5 
AIRS Version-5 minus CERES Aqua Edition 2.5 
CERES Terra minus CERES Aqua 


Figure 2a. Global monthly mean OLR differences among AIRS, CERES Terra, and 
CERES Aqua; and dashed line showing the average value of each OLR difference for 
the overlap periods between each data set. 
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September 2002 through June 2010 
Global Clear Sky OLR Time Series Differences (Watts/m 2 ) 



AIRS Version-5 minus CERES Terra Edition 2.5 
AIRS Version-5 minus CERES Aqua Edition 2.5 
CERES Terra minus CERES Aqua 


Figure 2b. Global monthly mean OLR C lr differences among AIRS, CERES Terra, and 
CERES Aqua; and dashed line showing the average value of each OLRclr difference 
for the overlap periods between each data set. 
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Outgoing Longwave Radiation (Watts/m 2 ) 
AIRS Version-6 minus AIRS Version-5 
7-Day Average 




GLOBAL ME AN= -7.064 



STANDARD DEV= 4.142 CORRELATION:: 0.994 


Figure 3a. The AIRS Science Team Version-6 OLR minus Version-5 OLR for the 1:30 
AM/PM average of the seven day period, September 6, 2002, January 25, 2003, 
September 29, 2004, August 5, 2005, February 24, 2007, August 10, 2007, and 
May 30, 2010. 
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Clear Sky Outgoing Longwave Radiation (Watts/m 2 ) 
AIRS Version-6 minus AIRS Version-5 
7-Day Average 




GLOBAL MEAN = -7.539 



STANDARD DEV= 3.845 CORRELAT I 0N= 0.993 


Figure 3b. The AIRS Science Team Version-6 OLR C lr minus Version-5 OLR C lr for the 
1 :30 AM/PM average of the same seven day period. 
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September 2002 through April 201 1 
Global OLR Anomaly Time Series 
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AIRS Version-5 

CERES Terra 

AIRS Version-5 minus CERES Terra 


Figure 4a. Global mean AIRS and CERES Terra monthly mean OLR anomaly time 
series and their anomaly differences for the overlap period September 2002 through 
June 2010. 
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September 2002 through April 201 1 
Tropical OLR Anomaly Time Series 



AIRS Version-5 
CERES Terra 

AIRS Version-5 minus CERES Terra 
AIRS El Nino Index multiplied by 3 


Figure 4b. Tropical (20°N to 20°S) mean monthly mean OLR anomaly time series and 
anomaly differences of AIRS and CERES Terra for the overlap period September 2002 
through June 2010. Also shown in black is the AIRS El Nino index multiplied by 3. 
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September 2002 through April 201 1 
Global Clear Sky OLR Anomaly Time Series 



AIRS Version-5 
CERES Terra 

AIRS Version-5 minus CERES Terra 


Figure 5a. Global mean monthly mean OLR C lr anomaly time series and anomaly 
differences of AIRS and CERES Terra for the overlap period September 2002 through 
June 2010. 
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September 2002 through April 201 1 
Tropical Clear Sky OLR Anomaly Time Series 



AIRS Version-5 
CERES Terra 

AIRS Version-5 minus CERES Terra 
AIRS El Nino Index multiplied by 2 


Figure 5b. Tropical (20°N to 20°S) mean monthly mean OLR C lr anomaly time series 
and anomaly differences of AIRS and CERES Terra for the overlap period September 
2002 through June 2010. Also shown in black is the AIRS El Nino index multiplied by 2. 
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September 2002 through June 2010 
OLR Zonal Mean Anomaly 
Average Rate of Change Time Series 



AIRS Version-5 
CERES Terra 

AIRS Version-5 minus CERES Terra 


Figure 6a. Zonal mean Average Rates of Change (ARC's) for AIRS and CERES Terra 
OLR monthly mean time series for the period September 2002 through June 2010, as 
well as the difference of the two sets of ARC’S. The majority of the decrease in global 
OLR during this period originates in the tropics south of 8°N. 
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September 2002 through June 2010 
Clear Sky OLR Zonal Mean Anomaly 
Average Rate of Change Time Series 



AIRS Version-5 
CERES Terra 

AIRS Version-5 minus CERES Terra 


Figure 6b. Zonal mean ARC'S for AIRS and CERES Terra OLR C lr monthly mean time 
series for the period September 2002 through June 2010, as well as the difference of 
the two sets of ARC’S. 
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AIRS Version-5 

OLR Anomaly Average Rate of Change (Watts/m 2 /yr) 
September 2002 through June 2010 


90 N 
60 N 
30 N 
EQ 
30 S 
60 S 
90 S 



-3.72 -3.00 -2.04 -1.32 -0.60 0.00 0.60 1.32 2.04 3.00 3.72 


GLOBAL ME AN= -0.088 STANDARD DEV= 0.597 


Figure 7a. Spatial 1° latitude by 1° longitude distribution of global AIRS OLR ARC’S over 
the time period September 2002 through June 2010. OLR Region 1 is outlined in gray 
and OLR Region 2 is outlined in black in this and most subsequent figures showing 
spatial distributions of ARC’S of different parameters. 
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CERES Edition-2.5 

OLR Anomaly Average Rate of Change (Watts/m 2 /yr) 
September 2002 through June 2010 



60 W 


0.00 0.60 1.32 2.04 3.00 


GLOBAL ME AN = -0.065 STANDARD DEV= 0.586 


Figure 7b. Spatial 1° latitude by 1° longitude distribution of global CERES OLR ARC’S 
over the time period September 2002 through June 2010. 
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AIRS Version-5 minus CERES Edition-2.5 
OLR Anomaly Average Rate of Change (Watts/m 2 /yr) 



n I i i I i r 

-4.35 -3.45 -2.55 -1.65 -0.75 0.00 0.75 1.65 2.55 


GLOBAL MEAN = -0.023 STANDARD DEV= 0.180 CORRELAT I 0N= 


3.45 4.35 

0 . 949 


Figure 7c. Spatial distribution of ARC’S of AIRS minus CERES OLR. 
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AIRS Version-5 

Clear Sky OLR Anomaly Average Rate of Change (Watts/m 2 /yr) 




-3.72 -3.00 -2.04 -1.32 -0.60 0.00 0.60 1.32 2.04 3.00 3.72 


GLOBAL ME AN= -0.004 STANDARD DEV= 0.214 


Figure 8a. Spatial 1° latitude by 1° longitude distribution of global AIRS OLRclr ARC’S 
over the time period September 2002 through June 2010. 
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CERES Edition-2.5 

Clear Sky OLR Anomaly Average Rate of Change (Watts/m 2 /yr) 
September 2002 through June 2010 
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EQ 


30 S 


60 S 


90 S 



-3.72 -3.00 -2.04 -1.32 -0.60 0.00 0.60 1.32 2.04 

GLOBAL ME AN= -0.069 STANDARD DEV= 0.233 


Figure 8b. Spatial 1° latitude by 1° longitude distribution of global CERES OLRclr 
ARC’S over the time period September 2002 through June 2010. 
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AIRS Version-5 minus CERES Edition-2.5 
Clear Sky OLR Anomaly Average Rate of Change (Watts/m 2 /yr) 




- 4.35 - 3.45 - 2.55 - 1.65 - 0.75 0.00 0.75 1.65 2.55 3.45 4.35 


GLOBAL ME AN= 0.065 STANDARD DEV= 0.187 CORRELAT I 0N= 0.762 


Figure 8c. Spatial distribution of ARC’S of AIRS minus CERES OLRclr. 
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AIRS OLR Anomaly (Watts/m 2 ) 

Tropics 5°N to 5°S 
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Figure 9a. Hovmoller diagram for time series of monthly mean AIRS OLR anomalies 
integrated over the latitude range 5°N through 5°S in each 1° longitude bin for the 
period September 2002 through April 201 1 . The longitudinal domain of OLR Region 1 
is indicated by the gray vertical lines in this and most subsequent Hovmoller diagrams. 
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CERES OLR Anomaly (Watts/m 2 ) 
Tropics 5°N to 5°S 

Monthlies, September 2002 throuah June 2010 
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Figure 9b. Hovmoller diagram for time series of monthly mean CERES OLR anomalies 
integrated over the latitude range 5°N through 5°S in each 1° longitude bin for the 
period September 2002 through June 2010. 
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AIRS minus CERES OLR Anomaly (Watts/m 2 ) 
Tropics 5°N to 5°S Correlation^. 993 
Monthlies, September 2002 throuah June 2010 
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Figure 9c. Hovmoller diagram for the difference between AIRS OLR and CERES OLR. 
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AIRS Clear Sky OLR Anomaly (Watts/m 2 ) 
Tropics 5°N to 5°S 

Monthlies, September 2002 through April 201 1 
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Figure 10a. Hovmoller diagram for time series of monthly mean AIRS OLR C lr 
anomalies integrated over the latitude range 5°N through 5°S in each 1° longitude bin 
for the period September 2002 through April 2011. 
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CERES Clear Sky OLR Anomaly (Watts/m 2 ) 
Tropics 5°N to 5°S 

Monthlies, September 2002 throuah June 2010 
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Figure 10b. Hovmoller diagram for time series of monthly mean CERES OLR C lr 
anomalies integrated over the latitude range 5°N through 5°S in each 1° longitude bin 
for the period September 2002 through June 2010. 
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AIRS minus CERES Clear Sky OLR Anomaly (Watts/m 2 ) 
Tropics 5°N to 5°S Correlation^. 91 3 
Monthlies, September 2002 through June 2010 
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Figure 10c. Hovmoller diagram for the difference between AIRS OLR and CERES 
OLRclr 
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AIRS Version-5 Surface Skin Temperature Anomaly 
Average Rate of Change (K/yr) 
September 2002 through April 201 1 
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Figure 11a. Spatial distribution of Average Rate of Change of AIRS surface skin 
temperature (K/yr) over the period September 2002 through April 201 1 . The spatial 
domain of the AIRS El Nino region is outlined in green. 
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AIRS Version-5 500mb Specific Humidity Anomaly 
Average Rate of Change (%/yr) 
September 2002 through April 201 1 
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Figure 1 1b. Spatial distribution of Average Rate of Change of AIRS 500 mb specific 
humidity (%/yr) over the period September 2002 through April 2011. 
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AIRS Version-5 OLR Anomaly 
Average Rate of Change (Watts/m 2 /yr) 
September 2002 through April 201 1 
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Figure 11c. Spatial distribution of Average Rate of Change of AIRS OLR (W/m 2 /yr) over 
the period September 2002 through April 201 1 . 
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AIRS Version-5 Clear Sky OLR Anomaly 
Average Rate of Change (Watts/m 2 /yr) 
September 2002 through April 201 1 
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Figure lid. Spatial distribution of Average Rate of Change of AIRS Clear Sky OLR 
(W/m 2 /yr) over the period September 2002 through April 2011. 
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AIRS Surface Skin Temperature Anomalies (K) 
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Figure 12a. Hovmoller diagram of AIRS Surface Skin Temperature anomalies (K). The 
longitudinal domain of the El Nino region is enclosed between the green vertical lines. 
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AIRS 500 mb Specific Humidity Anomalies (%) 
Tropics 5°N to 5°S 

Monthlies, September 2002 through April 201 1 
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Figure 12b. Hovmoller diagram of AIRS 500 mb specific humidity anomalies (%). The 
longitudinal domain of Region 1 is enclosed between the gray vertical lines. 
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AIRS Effective Cloud Fraction Anomalies (%) 
Tropics 5°N to 5°S 

Monthlies, September 2002 through April 201 1 
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Figure 12c. Hovmoller diagram of AIRS effective cloud fraction anomalies (%). The 
longitudinal domain of Region 1 is enclosed between the gray vertical lines. 
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AIRS Version-5 Regional Anomaly Time Series 
September 2002 through April 201 1 



OLR averaged over Region 1 lagged by 3 months 
OLRclr averaged over Region 1 lagged by 3 months 


Figure 13a. AIRS OLR and OLRclr anomaly time series averaged over OLR Region 1 
superimposed on the AIRS El Nino index multiplied by 4 and plotted 3 months earlier. 
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AIRS Version-5 Regional Anomaly Time Series 
September 2002 through April 201 1 



OLR averaged over Region 2 
OLRclr averaged over Region 2 


Figure 13b. AIRS OLR and OLRclr anomaly time series averaged over OLR Region 2 
superimposed on the AIRS El Nino index multiplied by 10 without a time lag. 
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